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Abstract: Wc propose a novel probabilistic method for detection of ob- 



jects in noisy images. The method uses results from percolation and random 
s ! graph theories. We present an algorithm that allows to detect objects of 

o 

I unknown shapes in the presence of random noise. The algorithm has linear 

complexity and exponential accuracy and is appropriate for real-time sys- 
^ tems. We prove results on consistency and algorithmic complexity of our 

procedure. 
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1. Introduction 

In this paper, we propose a new efficient technique for quick detection of objects 
in noisy images. Our approach uses mathematical percolation theory. 

Detection of objects in noisy images is the most basic problem of image analy- 
sis. Indeed, when one looks at a noisy image, the first question to ask is whether 
there is any object at all. This is also a primary question of interest in such 



diverse fields as, for example, cancer detection ( Ricci- Vitiani et al. (2007)), au- 



tomated urban analysis (Negri et al. (20061), detection of cracks in buried pipes 



(Sinha and Fieguth (2006)), and other possible applications in astronomy, elec- 
tron microscopy and neurology. Moreover, if there is just a random noise in the 
picture, it doesn't make sense to run computationally intensive procedures for 
image reconstruction for this particular picture. Surprisingly, the vast majority 
of image analysis methods, both in statistics and in engineering, skip this stage 
and start immediately with image reconstruction. 

The crucial difference of our method is that we do not impose any shape 
or smoothness assumptions on the boundary of the object. This permits the 
detection of nonsmooth, irregular or disconnected objects in noisy images, under 
very mild assumptions on the object's interior. This is especially suitable, for 
example, if one has to detect a highly irregular non-convex object in a noisy 
image. Although our detection procedure works for regular images as well, it is 
precisely the class of irregular images with unknown shape where our method 
can be very advantageous. 

Many modern methods of object detection, especially the ones that are used 
by practitioners in medical image analysis require to perform at least a prelim- 
inary reconstruction of the image in order for an object to be detected. This 
usually makes such methods difficult for a rigorous analysis of performance and 
for error control. Our approach is free from this drawback. Even though some 



papers work with a similar setup (see Arias-Castro et al. (2005)), both our ap- 



proach and our results differ substantially from this and other studies of the 
subject. We also do not use any wavelet-based techniques in the present paper. 
We view the object detection problem as a nonparametric hypothesis testing 
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problem within the class of discrete statistical inverse problems. 

In this paper, we propose an algorithmic solution for this nonparamctric hy- 
pothesis testing problem. We prove that our algorithm has linear complexity 
in terms of the number of pixels on the screen, and this procedure is not only 
asymptotically consistent, but on top of that has accuracy that grows exponen- 
tially with the "number of pixels" in the object of detection. The algorithm has 
a built-in data-driven stopping rule, so there is no need in human assistance to 
stop the algorithm at an appropriate step. 

In this paper, we assume that the original image is black-and-white and that 
the noisy image is grayscale. While our focusing on grayscale images could have 
been a serious limitation in case of image reconstruction, it essentially does not 
affect the scope of applications in the case of object detection. Indeed, in the 
vast majority of problems, an object that has to be detected either has (on the 
picture under analysis) a color that differs from the background colours (for 
example, in roads detection), or has the same colour but of a very different 
intensity, or at least an object has a relatively thick boundary that differs in 
colour from the background. Moreover, in practical applications one often has 
some prior information about colours of both the object of interest and of the 
background. When this is the case, the method of the present paper is applicable 
after simple rescaling of colour values. 

The paper is organized as follows. Our statistical model is described in details 
in Section [2j Suitable thresholding for noisy images is crucial in our method and 
is developed in Section [3] A new algorithm for object detection is presented in 
Section |4j Theorem [l] is the main result about consistency and computational 
complexity of our testing procedure. An example illustrating possible applica- 
tions of our method is given in Section [5j Appendix is devoted to the proof of 
the main theorem. 

2. Statistical model 

Suppose we have a two-dimensional image. For numerical or graphical processing 
of images on computers, the image always has to be discretized. This is achieved 

imsart-generic ver. 2007/04/13 file: Image_Analysis_and_Percolation_Square.tex date: January 13, 2013 



M. Langovoy and O. Wittich/ Detection in noisy images and site percolation 4 

via certain pixelization procedure. In our setup, we will be working with images 
that are already discrete. 

In the present paper we are interested in detection of objects that have a 
known colour. This colour has to be different from the colour of the background. 
Mathematically, this is equivalent to assuming that the true (non-noisy) images 
are black-and-white, where the object of interest is black and the background 
is white. 

In other words, we are free to assume that all the pixels that belong to the 
meaningful object within the digitalized image have the value 1 attached to 
them. We can call this value a black colour. Additionally, assume that the value 
is attached to those and only those pixels that do not belong to the object in 
the non-noisy image. If the number is attached to the pixel, we call this pixel 
white. 

In this paper we always assume that we observe a noisy image. The observed 
values on pixels could be different from and 1, so we will typically have a 
greyscale image in the beginning of our analysis. It is also assumed that on each 
pixel we have random noise that has the known distribution function F; the 
noise at each pixel is completely independent from noises on other pixels. 

Let us formulate the model more formally. We have an N x N array of 
observations, i.e. we observe N 2 real numbers {Yij}fj =1 - Denote the true value 
on the pixel 1 < i, j < N, by Irriij, and the corresponding noise by aEij. 

Therefore, by the above assumptions, 

Yij = Irriij + osij , (1) 
where 1 < i, j < N, a > 0, and 



1, if belongs to the object; 
Irriij = { (2) 
0, if does not belong to the object. 

To stress the dependence on the noise level a, we write our assumption on the 
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noise in the following way: 



Sij ~ F, Esij — 0, Vareij = 1 . (3) 

The noise here doesn't need to be smooth, symmetric or even continuous. More- 
over, all the results below are easily transferred to the even more general case 
when the noise has arbitrary but known distribution function F gen ; it is not 
necessary that the noise has mean and finite variance. The only adjustment 
to be made is to replace in all the statements quantities of the form P( ■ /a) 
by the quantities F gen ( • ) . The Algorithm 1 below and the main Theorem [I] are 
valid without any changes for a general noise distribution F gen satisfying ^ 
and 

Now we can proceed to preliminary quantitative estimates. If a pixel (i, j) is 
white in the original image, let us denote the corresponding probability distribu- 
tion of Yij by Po- For a black pixel we denote the corresponding distribution 
of by Pj. We are free to omit dependency of Po and Pi on i and j in our 
notation, since all the noises are independent and identically distributed. 

Lemma 1. Suppose pixel has white colour in the original image. Then for 
all y € R: 

P (Y ij >y) = l-F( 1 iY (4) 

where F is the distribution function of the standardized noise. 
Proof. (Lemma [l): By 

PoiYj >y) = l- P(oe l3 <y) = l-F(^ 

□ 

Lemma 2. Suppose pixel has black colour in the original image. Then for 
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all y € K: 




(5) 



Pi(Yij <y) = P(l + asij <y) = P(ae lJ <y-l) = F 



( 



y-i 



a 



) 



□ 



3. Thresholding and graphs of images 

Now we are ready to describe one of the main ingredients of our method: the 
thresholding. The idea of the thresholding is as follows: in the noisy grayscale 
image {ly}^ ,- =1 , we pick some pixels that look as if their real colour was black. 
Then we colour all those pixels black, irrespectively of the exact value of grey 
that was observed on them. We take into account the intensity of grey observed 
at those pixels only once, in the beginning of our procedures. The idea is to 
think that some pixel "seems to have a black colour" when it is not very likely 
to obtain the observed grey value when adding a "reasonable" noise to a white 
pixel. 

We colour white all the pixels that weren't coloured black at the previous 
step. At the end of this procedure, we would have a transformed vector of O's 
and l's, call it {Y We will be able to analyse this transformed picture 
by using certain results from the mathematical theory of percolation. This is 
the main goal of the present paper. But first we have to give more details about 
the thresholding procedure. 

Let us fix, for each N, a real number ccq(N) > 0, ccq(N) < 1, such that there 
exists 0(N) € K satisfying the following condition: 
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Lemma 3. Assume that |6|) is satisfied for some 9(N) € K. Then for the 
smallest possible 6(N) satisfying |6p it holds that 



=!-«„(«). (7) 

Proof. (Lemma [3]) : Obvious by Lemma [l] □ 

In this paper we will always pick ao(N) = cyq for all N E N, for some 
constant «o > 0. But we will need to have varying ao( •) for our future research. 
We are prepared to describe our thresholding principle formally. Let pf te 



be the critical probability for site percolation on Z 2 (see Grimmett ( 1999 1 for 
definitions). 

As a first step, we transform the observed noisy image {Yi,j}ij=i in the 
following way: for all 1 < i,j < N, 

1. If Yij > 0(N), set Yij :— 1 (i.e., in the transformed picture the corre- 
sponding pixel is coloured black). 

2. If < 0(N), set Y^ := (i.e., in the transformed picture the corre- 
sponding pixel is coloured white). 

Definition 1. The above transformation is called thresholding at the level 0(N). 
The resulting vector {Yi,j}ij = i of iV 2 values (0's and l's) is called a thresholded 
picture. 

Suppose for a moment that we are given the original black and white image 
without noise. One can think of pixels from the original picture as of vertices 
of a planar graph. Furthermore, let us colour these N 2 vertices with the same 
colours as the corresponding pixels of the original image. We obtain a graph G 
with iV 2 black or white vertices and (so far) no edges. 

We add edges to G in the following way. If any two black vertices are neigh- 
bours (i.e. the corresponding pixels have a common side), we connect these two 
vertices with a black edge. If any two white vertices are neighbours, we connect 
them with a white edge. We will not add any edges between non-neighbouring 
points, and we will not connect vertices of different colours to each other. 
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Finally, we see that it is possible to view our black and white pixelized picture 
as a collection of black and white " clusters" on the very specific planar graph 
(a square N x N subset of the Z 2 lattice). 

Definition 2. We call graph G the graph of the (pure) picture. 

This is a very special planar graph, so there are many efficient algorithms 
to work with black and white components of the graph. Potentially, they could 
be used to efficiently process the picture. However, the above representation of 
the image as a graph is lost when one considers noisy images: because of the 
presence of random noise, we get many gray pixels. So, the above construction 
doesn't make sense anymore. We overcome this obstacle with the help of the 
above thresholding procedure. 

We make 9(N)— thresholding of the noisy image {Y%,j}ij=i as in Definition 
[lj but with a very special value of 0(N). Our goal is to choose 9(N) (and 
corresponding a (N), see ^) such that: 



(9) 



where p s ^ te is the critical probability for site percolation on Z (see Grimmett 



( |1999fl , |Ke~sten| ( |1982| )). In case if both M and ([£) are satisfied, what do we get? 



After applying the 9(N)— thresholding on the noisy picture {Yijuj—n we 
obtained a (random) black-and-white image {Y Let Gat be the graph 
of this image, as in Definition [2j 

Since Gjv is random, we actually observe the so-called site percolation on 
black vertices within the subset of Z 2 . From this point, we can use results from 
percolation theory to predict formation of black and white clusters on Gn , as 
well as to estimate the number of clusters and their sizes and shapes. Relations 
^ and Q are crucial here. 

To explain this more formally, let us split the set of vertices V n of the graph 
Gat into to groups: V N = V}J n UV^ ut , where V^tlVff 1 * = 0, and V^ 1 consists of 
those and only those vertices that correspond to pixels belonging to the original 
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object, while V^ ut is left for the pixels from the background. Denote G™ the 
subgraph of Gat with vertex set V^ m , and denote G^ 1 ' the subgraph of Gat with 
vertex set V£ ut - 

If © and ^ are satisfied, we will observe a so-called supercritical percolation 
of black clusters on G™ , and a subcritical percolation of black clusters on G™' . 
Without going into much details on percolation theory (the necessary introduc- 



tion can be found in Grimmett ( 1999) or Kesten ( 1982 )), we mention that there 
will be a high probability of forming relatively large black clusters on G™, but 
there will be only little and scarce black clusters on G^'. The difference between 
the two regions will be striking, and this is the main component in our image 
analysis method. 

In this paper, mathematical percolation theory will be used to derive quanti- 
tative results on behaviour of clusters for both cases. We will apply those results 
to build efficient randomized algorithms that will be able to detect and estimate 
the object {Imij}^j =1 using the difference in percolation phases on G 1 ^ 1 and 

Gout 
N ■ 

If the noise level a is not too large, then ^ and (|9| are satisfied for some 
0(N) G (0,1). Indeed, one simply has to pick 6(N) close enough to 1. On the 
other hand, if a is relatively large, it may happen that ([8| and ^ cannot both 
be satisfied at the same time. 

Definition 3. In the framework defined by relations ([lJl-Q and assumptions 
(Al) - (A3), we say that the noise level a is small enough (or 1-small), if the 
system of inequalities ^ and Q is satisfied for some 0(N) G R, for all N G N. 

A very important practical issue is that of choosing an optimal threshold 
value 9. From a purely theoretical point of view, this is not a big issue: once 
^ and ([9]) holds for some 6*, it is guaranteed that after 6— thresholding we will 
observe qualitatively different behaviour of black and white clusters in or outside 
of the true object. We will make use of this in what follows. 

However, for practical computations, especially for moderate values of TV, the 
value of 6 is important. Since the goal is to make percolations on Vjy" and Vfi ut 
look as different as possible, one has to make the corresponding percolation 
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probabilities for black colour, namely, 



as different as possible both from each other and from the critical probability 
pf te . There can be several reasonable ways for choosing a suitable threshold. 
For example, we can propose to choose 9(N) as a maximizer of the following 
function: 

(i - - *») 2 + (i - f (*£?a=i) - n-) 2 , do) 

provided that (Jsj) and ([9| holds. Alternatively, we can propose to use a maximizer 
of 



4. Object detection 

We either observe a blank white screen with accidental noise or there is an 
actual object in the blurred picture. In this section, we propose an algorithm 
to make a decision on which of the two possibilities is true. This algorithm is 
a statistical testing procedure. It is designed to solve the question of testing 
Hq : Iij — for all 1 < i,j < N versus Hi : Iij = 1 for some i,j. 

Let us choose a(N) € (0, 1) - the probability of false detection of an object. 
More formally, a(N) is the maximal probability that the algorithm finishes its 
work with the decision that there was an object in the picture, while in fact 
there was just noise. In statistical terminology, a(N) is the probability of an 
error of the first kind. 

We allow a to depend on N; a(N) is connected with complexity (and ex- 
pected working time) of our randomized algorithm. 

Since in our method it is crucial to observe some kind of percolation in the 
picture (at least within the image), the image has to be "not too small" in order 
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to be detectable by the algorithm: one can't observe anything percolation-alike 
on just a few pixels. We will use percolation theory to determine how "large" 
precisely the object has to be in order to be detectable. Some size assumption 
has to be present in any detection problem: for example, it is hopeless to detect 
a single point object on a very large screen even in the case of a moderate noise. 

For an easy start, we make the following (way too strong) largeness assump- 
tions about the object of interest: 

(Dl) Assume that the object contains a black square with the side of size 
at least (pi m (N) pixels, where 



lim l ° Sa Z\ =0. (12) 



/T>o\ r Vim( N ) 

(D2) hm — — — = oo . 

JV->oo log TV 



(13) 

Furthermore, we assume the obvious consistency assumption 



Vim(N)<N. (14) 
Assumptions (Dl) and (D2) are sufficient conditions for our algorithm to work. 



They are way too strong for our purposes. It is possible to relax ( 13 ) and to 
replace a square in (Dl) by a triangle-shaped figure. 



Although conditions ( |12[ ) and (13 (are of asymptotic character, most of the 
estimates used in our method are valid for finite N as well. 

Now we are ready to formulate our Detection Algorithm. Fix the false detec- 
tion rate a(N) before running the algorithm. 

Algorithm 1 (Detection). 

• Step 0. Find an optimal 9(N). 
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Step 1. Perform 8(N)— thresholding of the noisy picture {Yij}^j =1 . 
Step 2. Until 

{{Black cluster of size <fi m (N) is found} 

or 

{all black clusters are found}}, 



Run depth-first search (Tarjan (1972)) on the graph Gn of 
the 6»(7V)-thresholded picture {Y itj }^ j=1 

• Step 3. If a black cluster of size fi m (N) was found, report that 

an object was detected 

• Step 4. If no black cluster was larger than ip im (N), report that 

there is no object. 

At Step 2 our algorithm finds and stores not only sizes of black clusters, but 
also coordinates of pixels constituting each cluster. We remind that 9(N) is 
defined as in Gjv and {Yij}^^ were defined in Section [3J and tp im (N) 



is any function satisfying (12). The depth-first search algorithm is a standard 
procedure used for searching connected components on graphs. This procedure 
is a deterministic algorithm. The detailed description and rigorous complexity 



analysis can be found in Tarjan ( 1972 1, or in the classic book Aho et al. ( 1975 ), 
Chapter 5. 

Let us prove that Algorithm 1 works, and determine its complexity. 

Theorem 1. Let a be 1-small. Suppose assumptions (Dl) and (D2) are satis- 
fied. Then 

1. Algorithm 1 finishes its work in 0(N 2 ) steps, i.e. is linear. 

2. If there was an object in the picture, Algorithm 1 detects it with probability 
at least (1 — exp(— Ci(a)ipi m (N))) . 

3. The probability of false detection doesn't exceed min{a(iV), exp(—C2(cr)ipi m (N))} 
for allN > N(a). 
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The constants C± > 0, Gi > and N(a) € N depend only on a. 

Remark 1. Dependence on a implicitly means dependence on 9{N) as well, but 
this doesn't spoil Theorem [I] Remember that we can consider 9(N) to be a 



function of a in view of our comments before ( 10 1 and ( 11 ). 

Theorem [l] means that Algorithm 1 is of quickest possible order: it is linear 
in the input size. It is difficult to think of an algorithm working quicker in this 
problem. Indeed, if the image is very small and located in an unknown place 
on the screen, or if there is no image at all, then any algorithm solving the 
detection problem will have to at least upload information about 0(N 2 ) pixels, 
i.e. under general assumptions of Theorem [T] any detection algorithm will have 
at least linear complexity. 

Another important point is that Algorithm 1 is not only consistent, but that 
it has exponential rate of accuracy. 

It is also interesting to remark here that, although it is assumed that the 
object of interest contains a ipi m (N) x ifi m (N) black square, one cannot use 
a very natural idea of simply considering sums of values on all squares of size 
<Pim{N) x (pi m (N) in order to detect an object. Neither some sort of thresholding 
can be avoided, in general. Indeed, although this simple idea works very well 
for normal noise, it cannot be used in case of an arbitrary, possibly irregular 
or heavy-tailed, noise. For example, for heavy-tailed noise, detection based on 
non-thrcsholded sums of values over subsquares will lead to a high probability 
of false detection. Whereas the method of the present paper can still work in 
many cases. 



5. Example 

In this section, we outline an example illustrating possible applications of our 
method. We start with a real greyscale picture of a neuron (see Fig. [l}. This 
neuron is an irregular object with unknown shape, and our method can be very 
advantageous in situations like this. 

Basing on this real picture, we perform the following simulation study. We 
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Fig 1 . A part of a real neuron. 



add Gaussian noise of level a — 1.8 independently to each pixel in the image, 
and then we run Algorithm 1 on this noisy picture. A typical version of a noisy 
picture with this relatively strong noise can be seen on Fig. [2j We run the 
algorithm on 1000 simulated pictures. As a result, the neuron was detected in 
98.7% of all cases. At the same time, the probability of false detection was shown 
to be below 5%. Now we describe our experiment in more details. 

The starting picture (see Fig. [I]) was 450 x 450 pixels. White pixels have value 
and black pixels have value 1. Some pixels were grey already in the original 
picture, but this doesn't spoil the detection procedure. As follows from Theorem 
[I] our testing procedure is consistent at least when ^ and ^ are satisfied, i.e. 
when 



l-$( - ] < p s c lte = 0.58..., (15) 



p s c < i-n — , (i6) 
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Fig 2. A noisy picture. 

where 9 — 0(450) is the chosen threshold and <J> is the distribution function of 
the standard normal distribution. In our case, we have chosen a default threshold 
— 0.5. As can be seen from considerations at the end of Section[3| this threshold 
is reasonable but not the most effective one. The thresholded version of Fig. [2] 
is shown on Fig. [3] 



With this choice of 0, the system of ( 15 ) and ( 16 ) amounts to 



Taking into account the symmetry of $, this is satisfied if and only if 



1 

2a 



$( — I > pf e = 0.58. 



As can be found from the last equation, our testing procedure is asymptotically 
consistent at least in the noise level range < a < 2. We have chosen a = 1.8 
in our simulation study. In practice, Algorithm 1 can be consistently used for 
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Pit 



i 



Fig 3. A thresholded picture. 



stronger noise levels, because in fact there is a numerically significant difference 
not only between subcritical and supercritical phases of percolation, but also 
within each of the phases. 

Suppose the null hypothesis is true, i.e. there is no signal in the original pic- 
ture. By running Algorithm 1 on empty pictures of size 450 x 450 with simulated 
noise of level a — 1.8 and 9 = 0.5, one can find that with probability more than 
95% there will be no black cluster of size 191 or more on the thresholded pic- 
ture. Due to an exponential decay of maximal cluster sizes, it is a safe bet to 
consider as significant only those clusters that have more than, say, 250 pixels. 
A different and much more efficient way of calculating <p(N) for moderate sizes 



of N is proposed in Langovoy and Wittich ( 2009 1 . 



For moderate sample sizes, the algorithm is applicable in many situations 
that are not covered by Theorem 1. The object doesn't have to contain a square 
of size 190 x 190 in order to be detectable. In particular, for noise level a = 1.8, 
even objects containing a 40 x 40 square are consistently detected. The neuron 
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on Fig. [I] passes this criterion, and Algorithm 1 detected the neuron 987 times 
out of 1000 runs. 
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Appendix. Proofs. 



This section is devoted to proofs of the above results. Some crucial estimates 
from percolation theory are also presented for the reader's convenience. 

Proof. ( Theorem [lj : 

Part I. First we prove the complexity result. Finding a suitable (approximate, 



within a predefined error) 9 from (10) or (11) takes a constant number of oper- 
ations. See, for example, Krylov et al. ( |1976 ). 

The 9(N)— thresholding gives us {Yij}^j =1 and Gjv in 0(N 2 ) operations. 
This finishes the analysis of Step 1. 



As for Step 2, it is known (see, for example, Aho et al. (1975), Chapter 5, 



or 



Tarjan (1972)) that the standard depth-first search finishes its work also in 
0(N 2 ) steps. It takes not more than 0(N 2 ) operations to save positions of all 
pixels in all clusters to the memory, since one has no more than N 2 positions 
and clusters. This completes analysis of Step 2 and shows that Algorithm 1 is 
linear in the size of input data. 
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Part II. Now we prove the bound on the probability of false detection. Denote 



p out (N) := l-F^iy (17) 

a probability of erroneously marking a white pixel outside of the image as black. 
Under assumptions of Theorem [l] p ou t(N) < p^ lte . 
We prove the following additional theorem: 

Theorem 2. Suppose that < p ut{N) < p^ te . There exists a constant C3 = 
C?,(jPout{N)) > such that 



P Pout (N)(F N (n)) < exp(-nC 3 { Pout {N))), for all n > <f im (N) . (18) 

Here [n) is the event that there is an erroneously marked black cluster of size 
greater or equal n, lying in the square of size N x N corresponding to the screen. 
(An erroneously marked black cluster is a black cluster on Gn such that each of 
the pixels in the cluster was wrongly coloured black after the 9 —thresholding.) 

Before proving this result, we state the following theorem about subcritical 
site percolation. 

Theorem 3. (Aizenman- Newman) Consider site percolation with probability p 
on Z 2 . There exists a constant \ s u e = ^site(po) > such that 



P P0 {\C\ > n) < e -" A «^(Po) t f or a u n > 1 



(19) 



Here C is the open cluster containing the origin. 
Proof. (Theorem |3| : See Bollobas and Riordan ( 2006 1 



□ 



To conclude Theorem [2] from Theorem [3j we will use the celebrated FKG 
inequality (see Fortuin et al.| ( 1971 1, or |Grimmett (1999), Theorem 2.4, p. 34; 
see also Grimmett's book for some explanation of the terminology). 
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Theorem 4. If A and B are both increasing (or both decreasing) events on the 
same measurable pair (Q,T), then P{A PI B) > P{A) P{B) . 



Proof. (Theorem [2|: Denote by C(i,j) the largest cluster in the N x N screen 
containing the pixel with coordinates (i, j), and by C(0) the largest black cluster 
on the N x N screen containing 0. By Theorem [3j for alH, j: 1 < i,j < N: 

P P o^(N)(\C(0)\>n) < e-^Wjw), ( 20 ) 



Obviously, it only helped to inequalities ( 19 ) and ( 20 ) that we have limited our 
clusters to only a finite subset instead of the whole lattice 1? . On a side note, 
there is no symmetry anymore between arbitrary points of the N x N finite 
square; luckily, this doesn't affect the present proof. 

Since { |C(0)| > n} and { \C(i, j)\ > n} are increasing events (on the mea- 
surable pair corresponding to the standard random- graph model on Gjv), we 
have that { |C(0)| < n} and {\C(i,j)\ < n} are decreasing events for all i, j. 
By FKG inequality for decreasing events, 



P Pou t (N)(\C(i, j)\ < n for all i,j,l< i,j<N) 



> 



II II P P^(N)(\C^j)\<n) > (by@) 

l<i,j<N 

> ( 1 _ e -nA.o.(p.««)) Jva 

We denote below by C% the "a out of 6" binomial coefficient. It follows that 
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21 



Pp out (N)(FN(n)) = P PoUN) (B(i,j), 1 < i,j < N : \C(i,j)\ > n) 
< i_ (i _ e ~nX s ^( Potlt)) N 2 



N 



fc=0 

N 2 



= E c 



k ~ 1 r k P~ n X site(Pout) k 



N' 1 



k=l 



N 2 e -n\, ite (p out ) + ^ N 1 e -n\ Bit e{Pout)^ 



because we assumed in (18) that n > ipi m (N), and logA^ = o(ifi m (N)). More- 
over, we see immediately that Theorem [2] follows now with some C3 such that 

< C 3 { Pout (N)) < \ stt e{Pout{N)). □ 

The exponential bound on the probability of false detection follows from 
Theorem [21 

Part III. It remains to prove the lower bound on the probability of true detection. 
First we prove the following theorem: 

Theorem 5. Consider site percolation on Z 2 lattice with percolation probability 
p > p s ^ te . Let A n be the event that there is an open path in the rectangle [0, n] x 
[0, n] joining some vertex on its left side to some vertex on its right side. Let M n 
be the maximal number of vertex-disjoint open left-right crossings of the rectangle 
[0,n] x [0,n]. Then there exist constants C4 = Ci(p) > 0, C5 = C$(p) > 0, 
Cq = Cq(p) > such that 

P p {A n ) >l~ne- c * n , (21) 



P p (M n <C 5 n) < e- C6 ™, (22) 
and both inequalities holds for all n > 1. 
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Proof. (Theorem [5]): One proves this by a slight modification of the correspond- 
ing result for bond percolation on the square lattice. See proof of Lemma 11.22 



and pp. 294-295 in Grimmett (19991. 



□ 

Now suppose that we have an object in the picture that satisfies assump- 
tions of Theorem [lj Consider any <pi m (N) x tpi m (N) square in this image. After 
9— thresholding of the picture by Algorithm 1, we observe on the selected square 
site percolation with probability 

P1 ,JN) : =i- F (^)> ! .f. 

Then, by (21 ) of Theorem[5j there exists C4 = Ci(pi m (N)) such that there will 
be at least one cluster of size not less than (fi m (N) (for example, one could 
take any of the existing left-right crossings as a part of such cluster), provided 
that N is bigger than certain N(pi m (N)) = N(a); and all that happens with 
probability at least 

\-ne~ Gin > l-e- C3n , 
for some C 3 : < C 3 < C 4 . Theorem [T] is proved. 

□ 
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